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Strongly interacting Fermi gases are of great 
current interest. Not only are fermions the most 
common particles in the universe, but they are 
also thought to have a universal thermodynamic 
behavior for strong interactions 0, 0, Efl. Re- 
cent experiments on ultra-cold Fermi gases pro- 
vide an unprecedented opportunityto test uni- 
versality in the laboratory @, H, M 0> 0|- In 
principle this allows - for example - the interior 
properties of hot, dense neutron stars to be in- 
vestigated on earth. Here we carry out a de- 
tailed test of this prediction. We analyze results 
from three ultra-cold fermion experiments involv- 
ing two completely distinct atomic species in dif- 
ferent kinds of atomic trap environments 0?0 ? S|- 
The data is compared with the pre dict ions of a re- 
cent strong interaction theory [9j, 1 1 Ot ] . Excellent 
agreement is obtained, with no adjustable param- 
eters. By extrapolating to zero temperature, we 
show that the experimental measurements yield a 
many-body parameter [3 ~ 0.59 ± .07, describing 
the universal energy of strongly interacting Fermi 
gases. 

Experiments on ultra-cold Fermi gases at micro-Kelvin 
temperatures are revolutionizing many areas of physics. 
Their exceptional simplicity allows tests of many-body 
theory in areas long thought to be inaccessible. The 
ability to widely tune the effective interaction between 
fermions via a broad Feshbach resonance in gases of 6 Li 
and 40 K has allowed resonance models proposed in high- 
Tc superconductivity theory 11 ill to be implemented with 
fermionic ultra-cold atoms [TJ|. Fortunately, the Pauli 
exclusion principle stabilizes the resultin g ex cited molec- 
ular states against collisional damping |13|. This has 
permitted the experimental observation of the smooth 
evolution of the Fermi gas from the attractive regime of 
Bardeen-Cooper-Schrieffer (BCS) superfluidity through 
to a regime where molecules form a Bose-Einstein con- 
densate (BEC) 0, [B, [H ELI E3, EH- On the cusp of 
this BCS-BEC crossover, there is a strongly interacting 
regime — the so-called unitarity limit which leaves 
the inter-atomic distance as the only relevant length 
scale. 

At this point, the gas is expected to exhibit a univer- 
sal thermodynamic behaviour, independent of any micro- 
scopic details of the underlying interactions. The ground 
state energy Eq of a homogeneous gas at zero tempera- 



ture should be proportional to the free Fermi energy, Ep. 
Thus, Eq = Ep(l + f3), where (3 is a universal many-body 
parameter. Substantial experimental efforts have been 
carried out to verify the existence of universality 0, 0|, 
though so far there has been no conclusive confirmation. 
This is mainly due to the lack of reliable thermometry in 
the strongly interacting regime [13] , leading to an uncer- 
tainty in the finite temperature corrections. 

The universal parameters estimated from energies at 
the lowest accessible (but unknown) temperature, range 
from = -.68 ± 0.1 to P = -.54 ± 0.05 @, i, Efl 
Tl\ . Since these are not at zero temperature, there is an 
unknown correction factor required to obtain the ground 
state energy, and hence the true value of (3. There is 
a similar range of estimated theoretical values, though 
a more precise value of (3 — —0.58 ± 0.01 was recently 
obtained from zero temperature quantum Monte Carlo 
simulations [201 ] . 

This situation has dramatically improved in the most 
recent thermodynamic measurements on strongly inter- 
acting Fermi gases of 40 K and 6 Li atoms 0, E|, which 
allow accurate estimates of the energy in the universal 
regime from the fermionic cloud size. In experiments on 
40 K carried out at JILA 0|, an adiabatic magnetic field 
sweep is used to compare measurements in the strongly 
interacting and weakly interacting regimes, so that the 
non-interacting temperature is also known from the cloud 
size after the sweep. An important conceptual advance 
of the Duke group [H| who use 6 Li, was the realization 
that this gives a model-independent measurement of the 
entropy in the strongly-interacting regime, thus allowing 
a precision test Q of theoretical predictions of universal 
thermodynamics. A different approach at Rice [ij], also 
with 6 Li, makes use of the detailed density distribution 
to estimate temperature and entropy. 

These ground-breaking investigations provide mea- 
surements accurate to the level of a few percent, which is 
an exceptional accuracy in this challenging field of ultra- 
low temperature physics. 

In this Letter, we re-analyze all the available experi- 
mental data from these three laboratories, thus obtain- 
ing the measured entropy-energy relation of two com- 
pletely different strongly interacting trapped Fermi gases. 
We compare this directly with a single universal the- 
oretical prediction. We use a diagrammatic approach 
based on functional path-integrals [3, |21| together with 



the local density approximation to treat the inhomoge- 



2 



3.0 
2.5 
~ h 2.0 
S 1.5 
5 1.0 h 
0.5 
0.0 




■ Duke Li expt. 

• JILA 40 K expt. 

O Rice 6 Li expt. 

Our theory 

Ideal gas 

J i I i 







S/(Nk B ) 



Figure 1: (Color online) Illustration of the universal 
thermodynamics of a strongly interacting Fermi gas. 

Comparison between theoretical predictions and experimen- 
tal measurements on the entropy-energy relation of strongly 
interacting Fermi gases in a harmonic trap. Inset highlights 
the low entropy region. The energy is in units of the Fermi 
energy, which is the highest single-particle energy level of a 
non-interacting gas in the same trap with the same number of 
fermions, TV. Brown solid curve is our theoretical prediction, 
and the black dash-dotted curve is the ideal gas energy. The 
experimental data is from Refs. [1,0] and [1], on 6 Li and 40 K 
fermionic atoms, with error-bars taken from the experimental 
papers using appropriate conversions. All results agree with 
a single, universal theoretical curve. 



neous trap. Below the superfiuid transition, our calcu- 
lations include pairing fluctuations, which are important 
in the strongly interacting regime, due to the onset of 
molecule formation. This approach is described in detail 
elsewhere 0|. Above threshold we use the well-known 
Nozieres-Schmidt-Rink (NSR) theory [H, 0, 0|. We 
demonstrate a quantitative test of these thermodynamic 
predictions by comparing experimental results on both 
types of atom with a single theoretical curve, as shown 
in Fig (Q]). There are no adjustable parameters, so this 
provides strong evidence for universality. 

Further, by using power-law extrapolation to estimate 
finite temperature corrections, we are able to show that 
all the experimental data give rise to a universal many- 
body coefficient of 8 ~ —0.59 ± 0.07. This is in excellent 
agreement both with the Monte Carlo results [20j, 8 = 
—0.58 ± 0.01, and our earlier theoretical prediction [9j| of 
8 « -0.599. 

We summarize the experimental procedures,as typi- 
fied by the JILA work using 40 K atomic gases [7J. Here, 
the strongly interacting gas is prepared in a harmonic 
trap at the Feshbach resonance, and the potential en- 
ergy is measured from the observed radius. Next, the 
magnetic field is swept adiabatically to a zero scatter- 
ing length field, and the potential energy at this field 



is again measured, which gives the non-interacting tem- 
perature. From this data, we obtain the entropy of the 
interacting gas 0, since the energy-entropy relation of 
a non-interacting Fermi gas is known. The total energy 
is also obtained, as it is twice the potential energy at 
the resonance, owing to the virial theorem 0. Fig. 2a 
gives the predicted potential energy ratio as a function 
of the non-interacting temperature in the presence of a 
harmonic trap, together with the experimental measure- 
ments. The observed reduction of the potential energy 
in the strongly interacting regime is theoretically repro- 
duced. Converting the non-interacting temperature and 
the potential energy into the total entropy and energy, 
respectively, we obtain the entropy-energy relation for 
40 K gas in Fig. 2b. We find an excellent agreement be- 
tween the experimental data and theoretical predictions 
below threshold. There is a small discrepancy just above 
the critical temperature of (T/Tp)° « 0.25 or the criti- 
cal entropy of S c ~ 2.2Nks, where we expect that the 
above threshold NSR theory may be less reliable. This 
effect is clearly visible in Fig. 2a, which gives the original 
experimental measurements. The nonlinear transforma- 
tion used to obtain the entropy-energy relation in Fig. 
2b means that conventional rectangular error-bars give 
only a qualitative indication of the uncertainties in this 
figure. 

Next, we discuss in greater detail the 6 Li experiment in 
Duke [8|], which has greater accuracy but involves some 
extra corrections due to the anharmonic trap used, and 
residual interactions in the weakly interacting cloud. The 
strongly interacting Fermi gas of N = 1.3(2) x 10 5 atoms 
is prepared in a shallow Gaussian (anharmonic) trap 
V(r) = U {1 - cxp[-m(uj 2 j_p 2 +lj^z 2 )/(2U )}} at a mag- 
netic field B = 840 G, slightly above the resonance po- 
sition Bo = 834 G. The coupling constant fe« = —30.0, 
where kp is the Fermi wave vector and a is the s-wave 
inter-atomic scattering length, is sufficient large to en- 
sure the onset of the universal thermodynamic behavior. 
Experimentally, the entropy of the gas is measured by an 
adiabatic passage to a weak interacting field B = 1200 G, 
where kpa = —0.75 and the entropy and temperature is 
known from the cloud size after the sweep. The energy E 
is determined model independently from the mean square 
radius of the strongly interacting fermion cloud (2 
measured at 840 G, according to the virial theorem 
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where Ep — (3-ZVw^Wz) 1 / 3 = ksTp is the Fermi energy 
for an ideal harmonically trapped gas at the trap center, 
and zp is defined by 3muj 2 zp = Ep. The correction 
factor 1 — k accounts for the anharmonicity in the shallow 
trapping potential Uo ~ 10Ep. 

Fig. 3a shows the bare experimental data on the ratio 
of the mean square axial cloud size at 1200 G to that at 
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Figure 2: (Color online) Analysis of the 40 K experiment 
in JILA. a, Experimental data on the potential energy E po t 
normalized to that in the non-interacting regime Ep 0t for 
a 40 K gas at unitarity, as a function of the non-interacting 
gas temperature (T/Tf)°. It converges to unity as expected 
(dashed green horizontal line). The experimental data are 
compared with our theoretical prediction for a harmonically 
trapped, strongly interacting Fermi gas. b, Comparison be- 
tween theory and experiment on the entropy dependence of 
the energy for a unitarity gas. The experimental entropy is 
calculated from the non-interacting temperature, with error- 
bars converted from the experimental data. The energy is 
double the potential energy E po t due to the virial theorem [5(. 



840 G, as a function of the energy at 840 G, as compared 
to the theoretical simulations with the same realistic pa- 
rameters, except that we use a resonance field B for the 
strongly interacting gas. Good agreement is found, with 
no free parameters. As before, there is a small discrep- 
ancy between the raw data and theoretical predictions, 
just above the critical energy. We have recalculated the 
entropy corrections due to residual interactions in the 
1200G cloud to improve the accuracy at the 1% level, by 
using an above-threshold NSR theory. Calibration of the 
entropy from the measured mean square axial cloud size 
at 1200 G using the theoretically predicted dependence 
of the entropy on the size (inset in Fig. 3b) leads to the 
comparison for the entropy-energy relation, as shown in 
Fig. 3b. The agreement is even more impressive. 
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Figure 3: (Color online) Analysis of the 6 Li experiment 

in Duke, a, Experimental data on the ratio of the mean 
square cloud size at 1200 G, (-z 2 ) 1200 to that at 840 G, (z 2 ) g40 
for a 6 Li gas, as a function of total energy, are compared to 
the theoretical simulations. The data is obtained by an adi- 
abatic passage from a strongly interacting field B — 840 G 
to a weakly interacting field B = 1200 G. The total energy is 
measured at 840 G by utilizing the virial theorem, and is nor- 
malized with respect to the Fermi energy of a non-interacting 
gas. The theoretical curve are calculated under the same pro- 
cedure and parameters, except that the starting field is at 
the Feshbach resonance B = 834 G. As shown by the dashed 
green horizontal line, the ratio converges to unity at high en- 
ergy as the gas becomes more ideal. The arrow points to the 
theoretically predicted energy at transition point, b, Exper- 
imental data on the entropy dependence of the energy for a 
strongly interacting gas at 840 G, compared to the predic- 
tions from our strong interaction theory. The result for an 
ideal gas is also plotted. The entropy and error-bars of the 
gas are calibrated from the experimental mean square axial 
cloud size (z 2 ) 1200 , using our theoretical dependence of the 
entropy on the cloud size (as shown in the inset), which should 
be extremely accurate in the weakly interacting regime. 



We can now describe the procedures used to obtain Fig 
(1), which illustrates the universal thermodynamic be- 
havior of a strongly interacting Fermi gas. Here we have 
plotted all the measured data in a single figure, and com- 
pared them with our prediction for the entropy depen- 
dence of the energy of a harmonically trapped, strongly 



4 



interacting Fermi g well as that of an ideal Fermi 

gas. The slight shift of the experimental data in the 
Duke experiment due to the anharmonicity of the trap 
has been corrected, by subtracting the (small) theoret- 
ical difference between a shallow Gaussian trap and a 
harmonic trap for the energy at the same entropy. We 
plot also a single data point from the 6 Li experiment 
in Rice for the energy at their lowest entropy. The 
agreement between theory and experiment is excellent 
for almost all the measured data. Exactly the same the- 
ory is used in all cases, with results from three different 
laboratories @, 0, Q . The universal thermodynamics of a 
strongly interacting Fermi gas is strikingly demonstrated, 
independently of which atomic species we compare with. 

Just above the critical entropy S c ~ 2.2Nks, for the 
superfluid-normal fluid phase transition, there is a sug- 
gestion of a discrepancy between theoretical predictions 
and these precise measurements. At this point the above- 
threshold NSR theory is least accurate [13] • The data 
here may even indicate a first order superfluid phase tran- 
sition. We note that the exact order of phase transition 
for a strongly interacting Fermi gas is still an open ques- 
tion [25], and merits further investigation. 

A key feature of current ultra-cold Fermi gas exper- 
iments is that the lowest attainable entropy is around 
S = 0.7NkB, which corresponds to a temperature of 
0.10 — 0.15Tf at unitarity. This nonzero entropy or tem- 
perature affects the precise determination of the univer- 
sal many-body parameter (3. To remove the temperature 
dependence, we assume that in the low entropy regime 
(below the phase transition), there is a power law de- 
pendence of the energy on the entropy: E — E Q oc S a , 
as anticipated theoretically. For non-interacting Fermi 
and Bose gases, the power law exponents are 2 and 
4/3, respectively. The thermodynamics of a unitary 
gas should lie between an ideal Fermi gas and an ideal 
Bose-Einstein condensate. The fitting procedure leads to 
E Q /(NE F ) = 0.48 ± 0.03 and E /(NE F ) = 0.48 ± 0.04, 
for the Duke and JILA experiments, respectively, with a 
similar power law exponent a = 1.7 ±0.4. Here the error 
bar accounts for the fitting error only. Using the rela- 
tion E Q /N = (3/4)(l + I3f/ 2 E F for a harmonic trap Q, 
this gives rise to (3 ~ —0.59 ± 0.07: which agrees fairly 
well with the most accurate quantum Monte Carlo sim- 
ulations [20j, (3 = —0.58 ± 0.01, and our theoretical pre- 
dictions [9], (3 « —0.599. Our theoretical power-law pre- 
diction is a — 1.5 in the LDA regime, which also agrees 
with experiment. 
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Methods. — We briefly explain our analytic the- 
ory. This is an approximate method using perturba- 
tion theory summed to all orders, since no exact re- 
sults are known. In the homogeneous gas case, it relies 
on the many-body T-matrix approximation to account 
for the effects of collective Bogoliubov- Anderson modes, 
and extends the standard NSR approach to the broken- 
symmetry state 0| . This amounts to considering the 
contributions of Gaussian fluctuations around the mean- 
field saddle point to the thermodynamic potential (with 
Nambu notation), 



5n = i^lndct 



where 

xu(Q) 



xn(Q) xi 2 (Q) 
xi 2 (Q) xn(-Q) 



Qu{Q - K)G n {K) 



(2) 



1 

2^' 



4irh 2 a 
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are respectively the diagonal and off-diagonal parts of the 
pair propagator. Here, Q = (q,iv n ), K — (k, iuj m ), and 
Y^k ~ (l an d k are wave vectors, v n and 

u) m bosonic and fermionic Matsubara frequencies, respec- 
tively), m is the fermion mass, T the temperature, and 
£k = fi. 2 k 2 /2m, Q\\ and Q12 are BCS Green's functions 
with a variational order parameter A. Together with the 
mean-field contribution 
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where the excitation energy E^ = [(et — /i) 2 + A 2 ] 1 / 2 and 
the Fermi distribution function f(x) — 1/(1 + e x / kE>T ), 
we obtain the full thermodynamic potential Vl = Slo 
All the observables are calculated straightforwardly fol- 
lowing the thermodynamic relations, once the chemical 
potential fi and the order parameter A are determined. 
Number conservation in the form n = —dVl/dfi is strictly 
satisfied, yielding an exact identity for the pressure P and 
energy density £ of a unitarity gas: P = 2/3£ (3j, la]- For 
simplicity, in our calculations we determine the order pa- 
rameter at the level of mean field, using the gap equation 
dflo/dA — 0. Part of our approach was also previously 



derived using a functional integral method [2lfl. In the 
case of the normal Fermi liquid with vanishing order pa- 
rameter, the usual NSR formalism is recovered [22]. 

This type of perturbation theory with bare BCS Green 
functions in the pair propagators constitutes the simplest 
description of the BCS-BEC crossover, including the es- 
sential pair fluctuations. More sophisticated approxima- 
tions with dressed Green functions in the pair propaga- 
tors, i.e., the GGq [26] and the fully self-consistent GG 
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schemes [27L l28|. have also been proposed. In Ref. [Io| . 
we performed a comparative study of these different ap- 
proximation schemes for a unitarity gas in the normal 
state. A related calculation in the superfluid phase has 
also been carried out recently [2^]. Compared with the 
latest path- integral Monte Carlo simulations (29I . l30| , our 
analytic perturbation scheme seems to be the optimal 
choice for the calculations of the type required for the 
entropy-energy relation. However, there is a small region 
around the critical temperature where none of the cur- 
rent calculations are reliable, and we see indications of 
this in the data. 

To include the effects of the trap, we employ the local 
density approximation by assuming that the system can 
be treated as locally uniform, with a position dependent 
local chemical potential /i (r) = // — V(r), where V(r) 
is the trapping potential. The local entropy and energy, 
calculated directly from the local thermodynamic poten- 
tial using thermodynamic relations, are then summed to 
give the total entropy and energy. We note that in the 
presence of a harmonic trap, the exact identity P = 2/3£ 
yields the virial theorem [5] , which states that the poten- 
tial energy of the gas is a half of its total energy. 
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